SOFT MATRIX MODELS AND CHERN-SIMONS PARTITION 

FUNCTIONS 



Abstract. We study the properties of matrix models with soft confining po- 
tentials. Their precise mathematical characterization is that their weight func- 
tion is not determined by its moments. We mainly rely on simple considera- 
tions based on orthogonal polynomials and the moment problem. In addition, 
some of these models are equivalent, by a simple mapping, to matrix mod- 
els that appear in Chern-Simons theory. The models can be solved with q 
deformed orthogonal polynomials (Stieltjes-Wigert polynomials), and the de- 
formation parameter turns out to be the usual q parameter in Chern-Simons 
theory. In this way, we give a matrix model computation of the Chern-Simons 
partition function on and show that there are infinitely many matrix models 
with this partition function. 
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Random matrix models ^1 |S] for a review) have attracted great amount 
of interest both due to its mathematical structure and for its manifold physical 
applications as well. When studying matrix models, one is especially interested in 
the solution of matrix models in the large N approximation, which corresponds to 
the planar Feynman diagrams E] • These models are characterized by the N by 
N matrix variable M and by an Hamiltonian H — TyV (M) . After diagonalization 
of the matrix, one is able to work in the eigenvalue space, and to consider the 
following expression for the partition function: 



with (3 — 1, 2 or 4, depending on the symmetry. Equivalently, from the perspective 
of random matrix theory one studies the joint probability distribution function for 
the N eigenvalues of the matrix ITH|. It has the well-known general form: 
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1. Introduction 
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Lo {x) and V{x) are named weight function and confining potential, respectively. If 
the elements of the random matrix are believed to be statistically independent from 
each other, one obtains the quadratic confinement potential V (x) = x^, leading to 
the Gaussian ensembles of random matrices 18 . In the usual physical applications 
of matrix models 8 , the confining potentials V {x) are such that the weight function 
uj {x) = exp[— y(a;)] is determined by the knowledge of all of its positive integer 
moments {7n}^o ' '^here 7„ = J a;"w [x) dx. In this paper, we want to study some 
general properties of matrix models with a weight function that does not satisfy 
this property. Interestingly enough, a very particular model in this category shows 
considerable interest in the context of Chern-Simons theory |23II14| and topological 

strings [2iiniiini. 

We briefly remind now the formalism of orthogonal polynomials. Recall that the 
relevant quantities, such as the density of states and the correlation functions, are 
obtained from the two-point kernel |18|: 

(1-3) KN{x,y)^^Lpk{x)ipk{y), 

k=0 

with ipk (x) — LO {x)pk {x), where pk {x) are the orthogonal polynomials associated 
to the weight function oj {x) .This is the fundamental object, due to the following 
classical result |18| : 

(1-4) Rk{xi,X2, Xk) = det [Kn {Xp, Xq)]g p^-^ 2,...,n ■ 

where the correlation functions are defined as: 

(1.5) i?fe(a;i,X2, ...,Xfc) EE — — - / dxk+i-: dx^P ixi,X2, .-^xn) ■ 

For Hermitian ensembles {P — 2), the density of states is the kernel at the origin: 



Af-l 

(1.6) p{x) = 5]^fe(x)2 

fe=0 

(1.7) = uj{x){pn{x)p'j^_j^{x)-p'j^{x)pn-i{x)), 

the second expression follows when applying the Christoffel-Darboux formula |^ 
I18| . For Gaussian ensembles for example, this quantity tends to the well-known 
semi-circle law in the limit N oo jl8| . In this work we shall consider Hermitian 
ensembles, and this expression for the density of states will be useful to understand 
relevant conceptual points in our discussion. 

As mentioned, we are concerned with the special properties of matrix models 
with an indeterminate weight function. Our discussion is mainly based on the 
classical moment problem (see ^H] for a recent review, that we will follow). The 
connection with random matrix theory can be readily guessed by the fact that only 
the moments play a role in the orthogonalization procedure. Therefore the two 
sets, the moments and the orthogonal polynomials, are essentially equivalent. 

The paper is organized as follows: In the next Section, we give a short introduc- 
tion to Chern-Simons theory and, focussing on these recent developments, we show 
how the matrix models in I15| are directly related with an Hermitian matrix 
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model with a log-normal weight function. In this way, one can do exact computa- 
tions through the associated orthogonal polynomials. To illustrate this, we com- 
pute the partition fimction of Chern-Sinions theory on with gauge group U{N). 
Interestingly enough, the use of g-deformed orthogonal polynomials readily leads 
to the natural parameter of Chern-Simons theory q = e^^-' = e™+'= . Therefore, 
the orthogonal polynomials approach is a method for non-perturbative solutions 
in Chern-Simons theory. In Section 3, we show that there are actually infinitely 
many matrix models with the same partition function. Finally, we show that this 
is a general feature of models with very weak confining potential. To conclude, we 
present a summary and possible directions for further research. 

2. Chern-Simons theory on S^: Matrix Model formulation 

We begin by outlining basic and well-known facts about Chern-Simons theory, 
emphasizing recent results and trends that connect Chern-Simons theory with ma- 
trix models ^ E]. As it is well-known, Chern-Simons theory is a topological 
quantum field theory whose action is built out of a Chern-Simons term involving 
as gauge field a gauge connection associated to a group G on a three- manifold AI. 
The action is: 

(2.1) S(A)^—f Tr ( ylAd.4 + -AaAaA I , 

47r7Af V 3 y 

with k an integer number. The natural associated observables are the correlators of 
Wilson loops, and its main interest come from the fact that these correlators lead 
to quantum-group polynomial invariants of knots and links |23| . For a review of 
the field since the seminal work see |14| . 

We mention now some few developments related to Chern-Simons theory. As 
advanced in |24| . Chern-Simons gauge theory has a string description in the sense 
of 't Hooft ■ Indeed, as it is well-known, gauge theories with the SU (N) group 
admit a large N expansion. In these expansions, correlators are expanded in powers 
of 1/N while keeping the 't Hooft coupling fixed t = xN, with x the coupling 
constant of the gauge theory. In the case of Chern-Simons theory, this large N 
expansion is reminiscent of a string theory expansion. This connection between 
Chern-Simons and topological strings was first pointed out by Witten [23 and has 
been extended in ^Jj . For example, in |25 it is shown that if one wraps N D-branes 
on M in T*M, then the associated topological A-model is a [/ (N) Chern-Simons 
theory on the three manifold M. 

On the other hand, regarding matrix models, in |10) it was already pointed out 
that the structure of the partition function of Chern-Simons theory on with 
gauge group SU{N), resembles the usual expression for the partition function of a 
one matrix model in terms of its associated orthogonal polynomials. Moreover, in 
[T5] , it is shown that the partition function of Chern-Simons theory on and with 
gauge group U{N), is given by: 



i<j ^ ' i=l 

that describes open topological A strings on T*S^ with TV branes wrapping (see 
the details in ^Hl)- These type of models have been further considered in In 
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it was shown that topological strings for i?-branes are equivalent to Hermitian 
matrix models, then the idea in [J is to obtain the results in ^Sj, by applying 
mirror symmetry to obtain _B-brane matrix models. We will comment further on 
the results in later on. Now, we mainly focus on the study of (|2.2|l . For recent 
reviews along these lines see [TBj . 

We remind that, as usual in Chern-Simons theory, the string coupling constant 
5s is related with the k in (|2.1|l by: 

27ri 

2.3) gs = 1 -■ 

^ ' ^ k + N 

As explained in 1151, the limit of the parameter gs — > 0, lead to the usual Gaussian 
Unitary ensemble. One can argue, for example, that in this limit, the Gaussian 
becomes a Dirac delta function, and therefore the results are independent of the 
correlation factor, or level repulsion. In [Tl I15|. this model is essentially studied 
through the consideration of averages (in a Gaussian Unitary ensemble) of the 
following quantities: 



(2.4) e-S'^=i'^'="^("), with flfc = and ak {u) = ^ (u, - u,, . 

This method leads to a great amount of perturbative information. Nevertheless, we 
note that the connection with the usual Hermitian matrix models -that is, of the 
type given by - is much simpler. Let us consider the following simple change 
of variables: 

(2.5) e"-=x„ 
then, (|2.2|l . reads: 

(2.6) 7 n(— .)^n-r^e-'°^^^-^^^^|^- 

i<j 1=1 

The factor Yli can be readily absorbed by considering a simple but 

remarkable property of the log- normal function u (x) = e~'°s^ ^i/'^gs . 

(2.7) uj {xq) = ^xuj {x) , 

Note that this functional equation is just a particular case of the following elemen- 
tary identity: 

(2.8) ^-fe^log^(c(-"^''')x) ^ ^-fc=(log=x+aV4fc^-^^Iog=.) 

Therefore, the Chern-Simons matrix model is directly related to an ordinary 
Hermitian ensemble with log-normal weight function lo {x) — e^'°s xi/2gs^ Equiv- 
alently, one could have used the change of variables exp (m^ + Ngs) = Xi. Notice 
that the term Ng^ is necessary to cancel, together with the term coming from the 
Jacobian, the contribution Yii^i = ni<j FIT- Actually, this method was al- 

ready employed in [5], where the same matrix model seems to appear in a very 
different context. One of our goals is to make contact between this model and 
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Chern-Simons theory. As we shall see, this will comes out rather easily, due to the 
existence of a closed system of orthogonal polynomials associated to the log-normal 
weight. These polynomials are q-deformed polynomials known as Stieltjes-Wigert 
polynomials ]|21,. This opens the possibility of a complete solution for the model 
(|2.2() of j^E] (see also (HI), using the classical results for Hermitian ensembles and 
these polynomials. Secondly, we shall show -giving explicit examples- that there 
exists infinitely many matrix models -deformations of 1)2. 6|1 - with exactly the same 
(Chern-Simons) partition function. 

2.1. Chern-Simons partition function through the matrix model compu- 
tation. Using the previous results and the Stieltjes-Wigert orthogonal polynomials 
pi ■ we can find a complete solution of the matrix model considered in |15|. and, in 
particular, for the partition function. For this, we can use the following well-known 
result for the partition function in terms of the orthogonal polynomials |H] : 
(2.9) 

N 



••• / n ^ (xO dx^ Y[ (x, - xif = j^J^ ^ = N\a^ 

'' i=l i<l lli=0 



i<l 

where the coefficient a, is 




N-i 



(2.10) p,{x)^a^x' + 

The first step is to compute the Z associated to the Stieltjes-Wigert orthogonal 
polynomials. The coefficients are given by j21| : 



(2.11) a, =g(^™^{(l-g)...(l-g^)} 

then 



-1/2 



(2.12) 

and ao = g^^^. Therefore: 



N-l 

(2.13) Z,^ = N\q-^/^\{q-^^i^~=){l~q')''-' 



N-l 



but we are interested in : 



(2.14) Z,i„h = / n $6"^ n (2 sinh 

= (2.)--e-^|n'^-^^"*^n( 



i—1 i<j 

Making the identification fe^ = ^ , and therefore q — e~^= . It is remarkable 
that we have naturally obtained, from this q deformed orthogonal polynomials, the 
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usual q parameter in Chern-Simons theory. The partition function with the usual 
Stieltjes-Wigert weight, corresponds to: 

N 



(2.15) Z,^ = l'f[dx,7T-'/^ke-''"'°s\-.)Y[{x,-Xjf 

r ^ 1 2( ) 



Therefore: 



2—1 i<j 



(2.16) Z.„,= (|l)"''^!e^--K-) n 

To make explicit contact with the typical expressions for the Chcrn-Simons partition 
functions (as in (10) for example) we make some transformations on the product 
term. 



N-l 



N-l N-l , .s N-] 



e 



N-l , N-j 

(2.17) = e^^s=^(^'-i)e^^(^-i). n (2sin^ 



The final expression for the partition function is |15j : 

-ig^Ar(7V=_i) 
(2.18) Z = - 



iV! 

Af/2 



N-l , .s N-j 



1 n 



and since 5s = -^ijq^ we finally find: 

JV-l , . s N-j 

(2.19) ^ = e^-"^(fc + iV)-"/^n 

It was already mentioned in [10) . that the structure of the partition function 
resembled very much the general expression for the partition function of a one ma- 
trix model in terms of the associated orthogonal polynomials. Here, based on the 
result in 15|, we are making this statement precise through the explicit computa- 
tion. Notice that the differences between the SU{N) and U{N) comes essentially 
from the contribution of the partition function of U{1) (see, for example, the first 
reference in for more details). 

3. Infinitely many matrix models with the same Chern-Simons 

partition function 

In this second part we shall show that the matrix models involved that appear in 
Chern-Simons theory present a novelty in comparison with most of known matrix 
models (like the ones that are relevant in 2D quantum gravity). Broadly speaking, 
its weight function is so broad that there is a non-uniqueness issue that leads to 
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very practical and concrete results for the matrix model: there are infinitely many 
models with identical properties. In particular, identical (Chern-Simons) partition 
function. We begin by briefly reviewing some introductory results in the theory of 
the moment problem. 

3.1. Relevant results from the moment problem. The two basic moment 
problems are 

3.1.1. Hamburger moment problem. Given a sequence of reals 70,71,... , when is 
there a measure, dp, on (—00, 00) so that 

/oo 
x"dp{x) 
-00 

and if such measure exists, is it unique ? 

3.1.2. Stieltjes moment problem. Given a sequence of reals 70, 7i, ... , when is there 
a measure, dp, on (0, 00) so that 

POO 

(3.2) 7„= / x^'dpix) 

JO 

and if such measure exists, is it unique ? 

In our case, we are not concerned with the question regarding existence, since 
we always refer to the weight function lo (x) and its associated moments. Then, 
the main point is to know whether we are in the determinate or indeterminate 
case. To fix ideas, let us mention that there is a natural boundary for functions 
of the form e"^*^^-*, with V{x) of polynomial form. In the Hamburger case, the 
boundary is given by V{x) = |a;|" with a = 1, while in the Stieltjes case it is given 
by V{x) = |a;|" with a = 1/2. 

Then, in the usual matrix models considered (as reviewed in 8 for example), 
one is essentially always in the determined case. 

Consider now the following generalization of the log-normal weight: 

(3.3) /e(a;) =e"'°s'="(l + ^9sin(27rlogx)), with i9 e [-1, 1] . 
The moments of this function are |2(Jj : 

/>oo 

(3.4) 7„=/ fe{x)x"dx = V^ei^""+^\ 

Jo 

Note that all the integer moments are completely independent of the parameter 
'd. This means that all the functions in the family /e(a;) have the same moments. 
Thus, they are all undetermined by them. Conversely, one can say that the set of 
moments 7„ = -^/TreJ^" is an indeterminate set. We present now a detailed 
consideration of this family of functions. The main point is the computation of its 
Mellin transform. 

First, note that for 6* = we have the log-normal distribution e^ Its Mellin 

transform is: 

(3.5) M[/(x);s]=e-'/4Erf(-0. 

The other part of the Stieltjes function is e~'°s ^ sin (27rlna;). From (|3.4(l . we 
know that all of its integer moments are zero. The most enlightening possibility is 
to compute its Mellin transform: 
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Figure 1 . The oscillatory term in fg {x) 



(3.6) / = / a;^-ie-'°s'Wsin(27rloga;)dx. 

Jo 

We use the change of variables y — — + logx. The integral becomes: 

/oc 
e-y"+i('"+') sin 
-oo 

= v^e3(^'+i)sin(7r(l + s)). 

The change of variables is the one employed by Stieltjes [50], with a generic, 
complex s instead of an integer k. Clearly, after the change of variables, since sin 
is a periodic and odd function, its value is zero at s = fc. Nevertheless, the result 
is even clearer when looking at the full Mellin Transform: the contribution of the 
additional function is zero only for integer values of s. 

As a consequence, for example, we have that any of the following functions have 
the same integer moments. 



27r U 



dy 




Figure 2. 6* = 1/2 
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Figure 3.0 = 1 (solid) and 6 = -1 (dashed) 



Consequently, it is manifest now that there are infinitely many matrix models 
with exactly the same partition function. Therefore, it turns out that the partition 
function of U (N) Chern-Simons theory on can be given by an infinitely many 
different matrix models. From what we have seen, one possibility is: 



N , / \ 2 



(3.8) ff[ e-"' /'^^ (1 + sin (27r {ujg, + N))) ^ H ^inh ■ 

1=1 i<j ^ 

as always, with Be [—1,1]. Notice that due to the correspondence between (|2.2|) 
and (|2.6|l through H2.5() , the density of states of (|3.8|l does not satisfy the property: 



(3.9) / pe {u)u"-du = jn, 

J —OQ 

(that is, the moments are different for each matrix model). But rather: 

/oo 
pg (e") e("+i)("+^s=)du = 7„, with neZ. 
-oo 

To extract further physical consequences of these results and its implications 
for the connection between matrix models and topological quantum field theories 
seems an interesting open question. 

4. On matrix models with soft potentials 

In this concluding Section we show that the properties exhibited by the Chern- 
Simons matrix model are typical of models with weakly confining -soft- potentials. 
The mathematically precise description is that the potential is such that e"^'^-' is 
a function that is not determined by its moments. This feature has already been 
discussed in the context of random matrix theory and also in the particular case of 
q-deformed matrix models ,17, . In this last case, the one we are mainly interested, 
the non- uniqueness is strong since the potentials are asymptotically V{x) ~ log^ x 
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for X oo (17j and we have already mentioned that we are in the soft regime when 
the confinement provided by the potential is weaker than V{x) ~ x. 

Employing basic results on indeterminate functions w (x) , and the classical ran- 
dom matrix formulae, one can prove that the density of states of a matrix model is 
an indeterminate function if its weight function is indeterminate by its moments. 
This can be readily shown if we consider the expression for the density of states (for 
any N) and Krein proposition [T^, that essentially says that if dx 
is convergent, then the corresponding moment problem is undetermined. An anal- 
ogous case holds for the Stieltjes case. 

From (11.6(1 we have: 

\ogp{x) = \oguj{x) + log {pn{x)p'j^^i{x) - p'j^ {x)pN-i{x)) 

(4.1) = \oguj{x) + log Qn{x). 

with Qn{x) a polynomial. We write the expression in the second way, to emphasize 
that wc just have the addition of a polynomial term. Then: 

(4.2) / — — da; < oo. 



and 

/x{l + x) 



(4.3) / }^^^dx<oo. 



Therefore, the convergence or divergence of f"" i dx and r„°° , , dx is ex- 

clusively given by the convergence or divergence of /^^ '-2S^dx and ^^i^dx, 
respectively. Then, taking into account Krein proposition, it turns out that the ma- 
trix models characterized by an indeterminate weight function, have different (but 
indeterminate) density of states but the same partition function. Thus, we have 
seen that the density of states is an indeterminate function. Let us show explicitly 
that their moments are all equal: 

(4.4) pe{x) =LJe{x)QN{x) , 
that is clearly 9 dependent and the moments: 

(4.5) Sn,9^ x"-LUe (x) Qn (x) dx ^ Sn- 



are 0— independent, since the polynomial Qn (x) is always the same and then we 
are lead to a sum of moments of the indeterminate function uJe (x) , each one of 
them 0— independent. Now, we can show that the partition functions are identical. 
We use the following expression: 



(4.6) Zn^9 = J pN{x)dx = J UJe (x) QN{x)dx = Zn- 

That is, the partition function as the 0-moment of the density of states. Alterna- 
tively, the usual expression for the partition function in terms of the orthogonal 
polynomials can be used |8): 
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where ai are the following coefhcients of the associated orthogonal polynomials: 
Pn {x) = unx^ + ... This expression is valid for any weight function and only 
depends on the associated orthogonal polynomials, that are the same for any 9. 

Then, for example, all matrix models with the weight functions H3.3|l . have the 
same orthogonal polynomials (essentially the Stieltjes-Wigert see above). Nev- 
ertheless, it is plain that each one will have a different density of states, since we 
have a different weight function (as we have seen in the previous figures). 

In any case, it is clear that approaching matrix models only from the perspec- 
tive of orthogonal polynomials, there is no a priori reason to choose the log-normal 
weight, since any member of the infinite family has identical orthogonal polyno- 
mials. Furthermore, the family of functions above discussed is by no means the 
only one. From our analysis with the Mellin transform, is clear that many other 
functions instead of sin, can lead to the same phenomena. In addition, examples 
such as |3j: 

2.-5/2 

(4.8) w(a;) = - -— — , < a; < oo, 

(-a;;?)^ {-q/x;q)^ 

can be relatively easily found in the literature. This weight function also leads to the 
Stieltjes-Wigert polynomials, that are actually terminating basic hypergeometric 
series and have intimate relationships with theta functions [H] . In this work, we have 
discovered their remarkable connection with Chcrn-Simons theory and quantum 
topological invariants. In general, it turns out that the set V of solutions to an 
indeterminate moment problem contains very different types of measures: measures 
fi € V with a C°°— density, discrete measures and measures which are continuous 
singular ^ . Explicit cases of these last types are found in 6 . 



5. Conclusions and outlook 



We have seen that the matrix models that appear in Chern-Simons theory are 
exactly solvable and that have many interesting properties, which are novel in 
comparison with the matrix models that appear in 2D quantum gravity. Therefore, 
the orthogonal polynomials technique constitute an interesting non-perturbative 
method. We have checked this for the simple case of 5''^ and gauge group U{N) and, 
interestingly enough, the use of g-deformed orthogonal polynomials readily leads 
to the usual parameter of Chern-Simons theory. It seems feasible, and certainly an 
interesting task, to try to extend this result to other geometries -like lens spaces- 
and to other gauge groups as well [T^ . 

On the other hand, the non-uniqueness property that we have established does 
not only seem to be interesting from a mathematical perspective, but it may have 
some physical meaning as well. We have emphasized the fact that the original 
Chern-Simons matrix model can be deformed in infinitely many ways, leaving the 
partition function always unmodified. Since this is a Chern-Simons partition func- 
tion, we know it is a topological invariant. A well-known result in topology says 
that the topological invariant of a three-manifold is independent of the triangula- 
tion employed to compute the invariant. Thus, one is tempted to try to interpret all 
these matrix models as representing different triangulations of the manifold. Thus, 
in a sense, it would be interesting to try to obtain an analogous interpretation to 
the role of more ordinary matrix models in 2D quantum gravity, but in three di- 
mensions. This is a bit speculative line of research, but it may still be worth. In 
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any case, the non-uniqueness properties are a distinctive mathematical feature of 
the models and the way the it appears seems also rather interesting in itself, since 
the resulting logarithmic oscillations are quite typical of physical models where a 
quantum group structure is present. A more detailed account of this will appear 
elsewhere. 
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